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Abstract 

Energy gap and wave function in thin films of topological insulator is studied, based on tight-binding model. It is 
revealed that thickness dependence of the magnitude of energy gap is composed of damping and oscillation. The 
damped behavior originates from the presence of gapless surface Dirac cone in the infinite thickness limit. On the 
other hand, the oscillatory behavior stems from electronic properties in the thin thickness limit. 
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1. Introduction 

Recently topological insulator has attracted much attention HI H SB Hi- Topological insulator has been firstly 
predicted in a graphene with spin-orbit interaction 1^, which consists of two copies of quantum Hall system, and 
shows a quantized spin Hall eff'ect if z-component of spin of electrons is conserved. Generally speaking, the present 
non-trivial system is characterized by Z2-index introduced by Fu, Kane, and Mele iHlsl], and has a gapless heli- 
cal edge mode where spin current protected by time-reversal symmetry flows spontaneously. However, quantum 
spin Hall (QSH) phase of graphene has not been observed experimentally since the spin-orbit interaction that is the 
driving force of topological insulator is much small. After that, HgTe/HgCdTe quantum well has been theoretically 
proposed as a candidate of two-dimensional topological insulator |@], and confirmed experimentally ifiol [ill [12I1. 
Topological insulators have been realized in three dimensional systems (HlQIlll], e.g., Bii-j^Sbjc alloy 1II6I1 . the bi- 
nary compounds Bi2Se3, Bi2Te3 HvlIislIlSl 20], and Tl-based ternary compound TlBiSes ||21,[22,[23|]. Furthermore, 
the quaternary compounds have also been theoretically predicted ||2^|25t]. All of these systems have a single helical 
Dirac cone on the surface. 

Nowadays, many exotic quantum phenomena are expected originating from surface states of three-dimensional 
topological insulators 



exotic qui 

mm 



27, 2; 



2911 . However, in the actual systems, sufl&cient amount of carriers remain in 
the bulk due to the difficulty of fabrication of samples ||3QIi . Then the system becomes metallic and it is difficult to 
classify physical properties specific to surface Dirac cone lH^, EH [s^] . To resolve this problem, several approaches, 
e.g., chemical doping and surface adsorption 113 311 have been performed. The another new app roach to control the 
carrier is to fabricate high quality thin films, where carrier control by gating is possible ISIlSD. But thin films may 
have the diff'erent electronic states from that of the bulk. Especially, the surface states have an energy gap due to the 
hybridization between the Dirac cones on top and bottom on the film induced by finite-size eff'ect. 

Based on above backgrounds, experimental studies of thin films of topological insulators have started [37l 
38ll . Besides this, there have been many theoretical studies based on continuous models il^liO,!!!^, first principle 
calculations IH, 111, IH, 111, li^, and tight-binding model calculation 1I42I1 . Although continuous model is simple, 
it is valid only for the long wavelength and low-energy limits. First principle calculation gives detailed electronic 
states of thin film. But it is difficult to analyze complicated phenomena, e.g., transport properties, disorder eff'ects, 
and quantum many-body problems. On the other hand, tight-binding approach is useful to calculate these interesting 
phenomena numerically, because many-body interaction and impurity eff'ects are easily taken into account. However, 
electronic properties of thin film of topological insulator have not been fully studied based on tight-binding model. 
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In the present paper, we study electronic properties of thin film of Bi2Se3 based on a tight-binding model focusing 
on the film-thickness dependencies of energy gap and surface states. It is revealed that the magnitude of energy 
gap is seriously influenced by material parameters. The paper is organized as follows. In section 2, we introduce a 
tight-binding model based on the Hamiltonian proposed by Refs. |[l7,[l8|]. In section 3, we calculate energy spectrum 
of thin film for various number of quintuple layers by changing material parameters. In section 4, we conclude our 
results. 



2. Model 



We use the eff'ective model derived in Refs. |ll7lll8|] as 
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with 



8(k) = Co + Cikl + C2(kl + k^y\ 
M(k) = Mo + Mik^^ + M2(kl + kb. 
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where k+ = kx ± iky, and the base is taken as (|+, t), I+, i), |-, T). i))^ in which ± and T (i) denote the parity 
eigenvalue and spin respectively. Let us introduce the lattice model only with nearest neighbor hoppings in a tetragonal 
lattice with substitution as 



kiGi sin kiGi, 
As a result, the bulk Hamiltonian is derived as 



(kai) 2(1 - cos kiGi). 
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where 



8(kx, ky) 
M{k) 

A+ikx, ky) 



Co + 2Ci(l - cos^^c) + 2C2(2 - cos^;^^^ - cos^3;a), (6) 
Mo + 2Mi(l - cos^^c) + 2M2(2 - cos k^a - coskyO), (7) 
Ao(sin kxa ± i sin kyo), (8) 

with Ui (a = ax = ay, c = a^) being the lattice constant along /(= x, z)-direction. The relation between the original 
parameters and those in the present model is 

Mo = Mo, Co = Co, Ml = Mi/c^ Ci = CJc^, 
M2 = M2/a^, C2 = C2/a^,Ao = Ao/a, Bo = Bole. (9) 

In the following, we express the present Hamiltonian in real space along z-direction perpendicular to the quintuple 
layers to focus on the surface states. Here, translational invariance is satisfied for the direction parallel to the quintuple 
layers i.e., x- and j-directions. Then kx and ky are good quantum numbers. We apply open boundary condition only 
along z-direction, i.e., the system is regarded as a one-dimensional chain for fixed {kx, ky). This condition corresponds 
to (1 11) cleavage surface of actual Bi2Se3 which is easily cleaved. The corresponding Hamiltonian is given as follows. 



H{kx, ky) = ^ elikx, ky)Ho(kx, ky)en{kx, ky) 

n=l 

N,-l 

+ ^ [el(kx, ky)Hien+i(kx, ky) + h.c] , 



(10) 



where A^^ denotes number of quintuple layers. It is noted that a lattice point n in the above Hamiltonian corresponds 
to position of a quintuple layer in the actual crystal structure. The on-site energy is given by 



Ho(kx, ky) = 
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with 



Sikx, ky) = Co + 2C\ + C2(2 - cos kxa - cos kyd), 
M(kx, ky) = Mo + 2Mi + M2(2 - cos kxa - cos kya), 



(12) 
(13) 



and the hopping between the nearest layers is as follows 
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Since the Hamiltonian H(kx, ky) has an inversion symmetry, it follows that [H, P] = 0, or equivalently PH(kx, ky)P = 



H(-kx, -ky), where the parity operator P is defined by 



A^7 



P = ^4(/:;„/:3;)diag[l,l,-l,-l]c^^+i_„(/:;„/:3;). 



(15) 



By using the parity operator P, we can derive the topological invariants v, which can be deduced from the parity of 
each pair of Kramers degenerate occupied energy band at the four time-reversal points at (Fi = (0, 0), T2 = (tt, 0), 
Fs = (0, tt), F4 = (tt, tt),) in the Brillouin zone, 



4 N, 



= nn<^2m(F.)iPi02m(Fj), 



(16) 
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where (pmi^a) is the eigenvector of the Hamiltonian //(F^,), and {(pmi^a) I P I (pmi^a)) (= ±1) is the eigenvalue of parity 
operator P. 



3. Results and discussions 

We numerically obtain the eigenvalues and eigenvectors of bulk and surface states, diagonalizing the Hamiltonian 
given by eq. (fTOb . The value of parameters Mo, M2, Aq, Co, C2 are the same as in Ref. ifisll with using a = 4.14A. The 
values of Mi and Ci are determined so that the eigen-energy at Z-point in Brillouin zone coincides with that of first 
principle calculation in Ref. ifisll . The value of Bq is chosen in order to fit the dispersion along F - Z line as well as 
possible. The indirect energy gap in the bulk Hamiltonian is located between -0.071 eV and |Mo + Col = 0.29 eV, as 
derived from eq. ©. 

Figure [T] shows the energy spectrum for a slab geometry in the cases of A^^ = 3, 5, 9, and 16. We can clearly see 
that eigenstates exist within the bulk energy gap. These states can be regarded as surface states, which we can directly 
confirm from its density distribution localized in the vicinity of surface, as shown in Figure [3j The surface states 
have a large magnitude of energy gap Eg ~ 0.033 eV for = 3 since the two wave functions localized at the top 
and bottom surfaces overlap significantly, (see (a) in Figure [T]). The magnitude of the present energy gap becomes 
small with the increase of A^^. ((b) and (c) in Fig. [TJ) For A^^ = 16, the resulting Eg is significantly reduced to be 
0.15 X 10"^ eV ((d) in Figure[T]). Moreover, it is noted that the shape of valence subband depends on A^^. For = 3 
((a) in Fig. [T]), there are two valence subbands in -0.4 eV < £^ < OeV. The upper subband consists mainly of surface 
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Figure 1: Energy spectra of the bulk and surface states near F-point in the slab geometry for A/^ = 3, 5, 9, and 16 quintuple layers. The material 
parameters are set as Mq = -0.28 eV, Mi = 0.216 eV, M2 = 2.60 eV, Aq = 0.80 eV, Bq = 0.32 eV, Co = -0.0083 eV, Ci = 0.024 eV, C2 = 
1.77 eV, a = 4.14 A, and c = 9.55 A. 
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Figure 3: Density distribution p(z) of surface state at (kxa,kya) = (;r/32,;r/32) for = 16. There are two surface states located at the top and 
bottom in the system, z denotes the position of layer. The guided lines are drawn for view-ability. The material parameters are the same as in 
Figure [T] 



States since it is located in the bulk energy gap. The lower one, that is bulk energy band, is located at 0.3eV. The 

new subband appears between these two subbands at 0.2eV for = 5 ((b) in Fig. [T]). Simultaneously, valence 

subbands have a local minimum at ^ = 0, and an indirect energy gap is generated. For = 9 ((c) in Fig. [T]) there 
are much more subbands. The energy bands for A^^ = 16 as shown in (d) in Fig. [T]is almost similar to that of bulk 
three-dimensional topological insulator with ^ oo. 

In Figure [21 A^^ dependence of Eg is plotted. Hamiltonian of the monolayer system with A^^ = 1 is given by 
Ho(kx, ky) (see eq. (fTTl)). which is equivalent to that of HgTe/HgCdTe quantum well |9], and Eg is given by 2|Mo + 
2Mi\{= 0.30eV). Eg for A^, = 2 is also derived analytically as |2Mi - [5^ + 4(Mo + 2MifYI\= 0.0094eV). For 
A^^ = 2, 3, and 4, Eg decreases roughly exponentially as a function of A^^, and becomes Eg ~ 10"^ eV. Eg becomes 
much smaller than room temperature for A^^ > 5. It is also noted that Eg has an oscillatory behavior as a function A^^ 
whose period is almost 3. The similar behavior has been obtained based on continuous models 139,^, 41J and first 

principle calculations Elilliilillii]. 



Next, we investigate the relation between the magnitude of energy gap and wave functions. Figure [3] shows the 
density distribution of surface states for A^^ = 16, which is defined by 

P{Z) = {cl(k,,ky)C,(k,,ky)) , Z = 1, • • • , A^„ (17) 

where the expectation value is evaluated for the surface state with momentum {kxa.kyo) = (tt/ 32, tt/ 32). The solid 
(dashed) line denotes the density distribution of surface state located on the top z = I (bottom z = 16). The density 
distribution decays exponentially with oscillation whose period is nearly 3 quintuple layers. This period is almost the 
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Figure 4: Phase diagrams of topological insulator for bulk limit = oo and for monolayer = I. STI (WTI) denotes strong (weak) topological 
insulator for = oo. QSH and OI denotes quantum spin Hall insulator where spin Hall conductance is quantized and ordinary insulator respectively 
for A^^ = I. Ml = 0.016,0.070,0.116,0.140,0.216,0.316, and 0.41 6eV for (a), (b), (c), (d), (e), (f), and (g), respectively. These values correspond 
to those used in Figure [5] 




Figure 5: The magnitude of energy gap Eg as a function of for different values of M\. The closed (open) circle describes a two-dimensional 
topological invariant v = 0(y = 1). At the case (d). Eg vanishes for odd numbers of layers: = 1, 3, • • • , 15. 



same as that of A^^ dependence of Eg. Since the two wave functions located at the top and bottom surfaces oscillate 
spatially, the resulting Eg due to overlap between them also oscillates as a function of A^^. 

In the following, we focus on the material parameters dependencies of Eg. For simplicity, we fix all parameters 
except for Mi. Here, we choose seven cases of Mi as shown in Figured In order to understand electronic properties 
for the corresponding seven cases we have chosen, we show the phase diagram of the system for A^^ = 1 and A^^ = oo in 
Figure|4j In the limit for A^^ = cx), the system becomes weak topological insulator (WTI) for Mi < -Mo/4 = 0.07eV 
while it becomes strong topological insulator (STI) for Mi > -Mo/4 = 0.07eV. On the other hand, in the limit for 
A^^ = 1 the present system is QSH for Mi < -Mo/2 = 0.14eV, while ordinary insulator (01) for Mi < -Mo/2 = 
0.14eV. 

Figure [5] shows Eg for various values of Mi . The curve in Figure [T] coincides with the curve (e) obtained for 
Ml = 0.216eV in Figure [S] which is the same parameter as that of Bi2Se3. The curve (e) has a three-fold periodic 
damped oscillation. For Mi = 0.016eV (case (a) in Fig. O, where the system is WTI (QSH) with = oo (Nz = 1), 
Eg does not decay with the increase of A^^ since there is no gapless surface Dirac cone at f -point for = oo. For 
Ml = O.OTOeV (case (b) in Fig. O, which is the transition point between WTI and STI, where closing of the bulk 
energy gap occurs at F-point, Eg decreases monotonically as a function of A^^. For Mi = 0.1 16eV (case (c) in Fig. O, 
Eg decays exponentially except for A/^ < 5 as a function of A^^, since surface Dirac cone is generated with Nz = oo. 
Eg has a strong oscillation for Mi = 0.140eV (case (d) in Fig. O, where transition between 01 and QSH occurs for 
A^^ = 1. Eg becomes exactly zero for odd numbers of A^^ (See Appendix). The damped oscillation with four-fold 
periodicity appears at Mi = 0.316eV (case (f) in Fig. [5]) and at Mi = 0.416eV (case (g) in Fig. O. When the system 
is QSH for A^^ = I, Eg decreases monotonically in the wide parameter range of A^^. On the other hand. Eg shows 
a damped oscillation as a function of A^^ for Mi > 0.140eV, i.e., the system is 01 in the thin thickness limit. As 
we have seen above, the period of oscillation depends on Mi . It can be concluded that the A^^ dependence of Eg is 
sensitive to the material parameter Mi . We also show the topological invariant v in FigJSl The closed circle expresses 
non-topological phase with y = 0, while the open circle expresses topological phase with y = 1. For cases with (a), (b) 
and (c), non-topological phase emerges for even number of layers and topological phase emerges for odd number of 
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layers. Topological phase and non-topological phase appear oscillatory also for cases with (e), (f) and (g). However, 
the period of oscillation becomes three or four. These results are consistent with those by Liu et al. 11421] . Since the 
magnitude of energy gap becomes zero at the boundary between topological and non-topological phases, the period 
of oscillation of topological number y coincides with that of the energy gap. 



4. Conclusion 



We obtain the energy spectrum of surface states in a topological insulator based on tight-binding model. It is 
clarified that there are various types of thickness dependencies of E^. The origin of the dumped oscillatory behavior 
of is partitioned into two parts. The dumped behavior appears when the gapless surface Dirac cone is realized in 
the limit of A^^ = oo. The oscillatory behavior of E^ becomes prominent when the system approaches to 01 regime for 

Based on these results, we would expect various types of thin films by controlling material parameters, which 
could be controlled by external pressure along z-direction (c-axis) for Bi2Se3. Tuning of material parameters may be 
much more easier for optical lattices made from cold atoms 147|, Hsl [13, EH 113] • If we can tune the corresponding 
material parameters with the case (d) in Fig.4, the strong even-odd eff'ect is expected. In this case, various transport 
properties are sensitive to external fields. There are several future unresolved problems. The Anderson localizations 
of thin films have been recently studied |,38|, i53i [SJ] from various aspects. It is interesting to study this problem 
with various types of thin films with diff'erent electronic properties of topological insulator with diff'erent material 
parameters. 
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Appendix 

As shown in Fig. [51 the energy gap is exactly zero for the odd number of layers, and non-zero for the even number 
of layers at the case (d) in Fig. IH In the following, we derive this behavior of the energy gap. The Hamiltonian for A^^ 
at f point reads 



with 4A^^ X matrices E^^ and G^^ being 



^2 



= Yj 4diag [(Co + 2Ci), • • • , (Co + 2Ci)] q 

N,-l 



(18) 
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(20) 



where Mq -h 2Mi = at the case (d). The eigenvalue of H^X^^ 0) equals to Co + 2Ci -h E^ where E^ is the eigenvalue 
of Gn^ . If matrix G^^ has zero eigenvalue, the energy gap closes because G^^ has particle hole symmetry . We show 
\GnJ = detGA^^ = for odd number of A^^ at f point, as follows. 
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Thus the energy gap at the case (d) is exactly zero for the odd number of layers, and non-zero for the even number of 
layers, and oscillates strongly as a function of N^. 
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